Strain-Induced Conduction Band Spin Splitting in GaAs from First Principles 

Calculations 
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We use a recently developed self-consistent GW approximation to present first principles calcu- 
lations of the conduction band spin splitting in GaAs under [f fO] strain. The spin orbit interaction 
is taken into account as a perturbation to the scalar relativistic hamiltonian. These are the first 
calculations of conduction band spin splitting under deformation based on a quasiparticle approach; 
and because the self-consistent GW scheme accurately reproduces the relevant band parameters, it 
is expected to be a reliable predictor of spin splittings. We also discuss the spin relaxation time 
under [lfO] strain and show that it exhibits an in-plane anisotropy, which can be exploited to obtain 
the magnitude and sign of the conduction band spin splitting experimentally. 

PACS numbers: 71.70.-d, 71.70.Ej, 71.15.-m, 71.15.Qe ,71.15.Mb ) 



I. INTRODUCTION 

The increasing prospect of utilizing spin electronics 
with conventional semiconductors, calls for quantitative 
predictions of the spin relaxation of electrons in these 
materials*^ In semiconductors without inversion sym- 
metry, the spin relaxation rate is related to the relativis- 
tic splitting in the conduction band, an effect which also 
induces spin precession and is relevant for spin trans- 
port and inject ion In zinc-blende semiconductors, 
which are the most promising for spintronic applications, 
it is widely accepted that D'yakonov-Perel' (DP)1 is the 
dominant spin relaxation mechanism. Generally speak- 
ing, this mechanism is present when the spin degeneracy 
of the conduction band is lifted. The spin splitting can 
be viewed as a k-dependent effective magnetic field fi(k) 
which under certain scattering conditions relaxes the av- 
erage spin of the ensemble. The strength of the effective 
field depends on the material. In the general form we can 
add this field to the Hamiltonian as an effective Zeeman 
term, H(k) = l/2er • fi, where fl(k) (with the dimen- 
sions of energy) is proportional to B e ff(k). In the zinc- 
blende crystal structure there is no inversion symmetry; 
this leads to an l~i field with components^ 



n l D = 2~/ki(k. 
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where i = x, i + 1 = y,i + 2 = z, the indices obey a 
cyclic relationship (i + 3 = 1), and 7 is a constant that 
depends on the bulk properties of the material. This ef- 
fective field was first introduced by Dresselhaus^. In uni- 
axially deformed crystals there is an additional effective 



^stress — C{ti,i+lh+l ~ £1,1+2^4+2) + 
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ey is the strain tensor, C and B are material-dependent 
constants. The first part of the effective field originates 
from off-diagonal components of the strain tensor while 
the second part originates from diagonal ones. The sec- 
ond part appears only due to the spin orbit mixing of p 
and d states and therefore should be much weaker than 
the first part 1 ^, but a numerical estimation of B has not 
yet been performed, probably because of uncertainties 
concerning the T\2 states. 

The data for the values of 7, C and B in various zinc- 
blende semiconductors are very sparse. The most studied 
case is GaAs. However, the experimental data for 7 show 
a wide range of values, between 11.0 and 34.5 eV-A 3 ^ 
Theoretical calculations of this parameter also show a 
wide range of predicted values. Calculations based on 
kp method predict a value between 25 and 30 eV-A 3 * 1 ^ 
However, a first principles calculation by Cardona et al^ 
predicted a value of 15 eV-A 3 . Our recent first principles 
calculation^ predicted a value of 8.5 eV- A 3 , a lot smaller 
than the commonly cited value of 27.5 eV-A 3 . Recently, 
Krich et al^ used a semi-classical approach to estimate 
the effect of flu on the mean and variance of the conduc- 
tance in closed quantum dots and compared the results 
of their model with the experiment in Ref. |19| . They were 
able to reach a good agreement only when they used our 
value of 7, suggesting that the value of 7 in GaAs must 
be around 9 eV-A 3 . 

Also our current knowledge of C and B in GaAs is far 
from being satisfactory. The magnitude of C was esti- 
mated experimentally in Refs. [2^ and [2l| to be 8.1±2.5 
eV-A and 3.9 eV-A, respectively. The values of C cal- 
culated with Linear Combination of Atomic Orbitals 
(LCAO) and pseudopotentials are 3.75 and 11.24 eV-A, 
respectively^. The calculations are able to define also 
the sign of C; in Ref. [Tjj] it was found to be opposite to 
that of 7i22i No other attempt has been made since to 
calculate the sign of C. To our knowledge the magnitude 
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FIG. 1: The left panel represents the applied deformation for 
e < and the right for e > . 

and sign of B has not been estimated experimentally nor 
theoretically so far. Although, M. I. D'yakonov et al^- 
showed that the spin relaxation time has a very weak de- 
pendence on applied pressure upon application of [100] 
strain, possibly indicating that the magnitude of B is 
negligibly small. 

In this work we use a recently developed ah inito 
method based on the GW approximation to predict these 
parameters for GaAs. We will try to answer the impor- 
tant questions about the strengths and signs of the spin 
splittings caused by the two mechanisms, Eqs. Q] and [21 



II. METHOD 

The GW approximation can be viewed as the first term 
in the expansion of the non-local energy dependent self- 
energy S(r, r, uj) in the screened Coulomb interaction W. 
From a more physical point of view it can be interpreted 
as a dynamically screened Hartree-Fock approximation 
plus a Coulomb hole contribution^ 3 .. It is also a prescrip- 
tion for mapping the non-interacting Green function to 
the dressed one, G° — > G. In the Quasiparticle Self- 
consistent GW (QSGW) method a prescription is given 
on how to map G to a new non-interacting Green func- 
tion G — > G°. This is used for the input to a new itera- 
tion; we repeat the procedure G° — > G — > G° — > .... until 
converegence is reached. Thus QSGIL is a self-consistent 
perturbation theory, where the self-consistency condition 
is constructed to minimize the size of the perturbation. 
QSGW is parameter-free, independent of basis set and 
of the LDA2i. The method is described in great detail 
in Refs. HJ and [H. It has been shown that QSGW re- 
liably describes the band structure in a wide range of 
materialaiZiSSiSLSS. 

The QSGW method in the current implementation 
uses the Full Potential Linear Muffin Tin Orbital (FP- 
LMTO) metho d 29 ' 30 , so we make no approximations for 
the shape of the crystal potential. The smoothed LMTO 
basis includes orbitals with I < l max = 5 and both 3d and 
Ad are included in the basis. Ad are added in the form 
of local orbitals 25 -an orbital strictly confined to the aug- 
mentation sphere, which has no envelope function at all. 
As QSGW gives the self-consistent solution at the scalar 



TABLE I: Important band parameters for GaAs. Eo = 
E(F C 6 ) - E(T%) and E' = E(T$) - E(Tl) are the energies of 
the first two conduction bands at the L-point. Aso and Ag Q 
are the spin-orbit splittings between Tg and L7 for valence 
and conduction bands, respectively, rn^/m is the conduction 
band effective mass at L. Energies are in eV; 7 is in eV-A 3 , 
C and B are in eVA. An asterisk in front of the value indi- 
cates that this is a calculated value from another theoretical 
method. 
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"From Ref. |32] 
"From Ref. Ill 
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Trom Table IX of Ref. |l6t where only the absolute value is re- 
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■'From Ref. [LH ; where only the absolute value is reported 
9 Calculated with Pseudopotentials; From Table IX of Ref. [l6| 
^Calculated with LCAO; From Table IX of RefJlJ 
'Calculated with Pseudopotentials; From Ref. [LH 
^Calculated with LCAO; From Ref.QJ 

^Calculated with the three-band kp method; From Ref. I 111 

relativistic level, we add the spin-orbit operator, Hso, as 
a perturbation (it is not included in the self-consistency 
cycle). 

It has also been shown that the QSGW method sys- 
tematically overestimates the fundamental band gap in 
semiconductors by an amount of a few tenths of an eV, 
independent of the magnitude of the gap2S. This error 
is related to the fact that the vertex correction is not 
taken into account in the method and when taken into 
account a nearly perfect agreement with experiment is 
achieved 3 ^. Here, in order to obtain highly accurate re- 
sults with less computational effort, we take a simple 
but somewhat heuristic approach to correct the error. 
We considered a 'hybridized' QSGPL+LDA Hamiltonian 
with 

H a + H so - #lda + (1 - a)(E - V£ DA ) + H so . (3) 

In Ref. [l7| we found that for all IILV and II-VI semi- 
conductors studied a value of a w 0.2 gives excellent 
agreement of calculated band gap and other important 
band parameters with experiment. 

All band parameters presented in Table I, except C 
and B, are calculated for the undistorted lattice struc- 
ture. All signs are presented with the convention that the 
anion is at the origin and cation at (0.25,0.25,0.25). We 
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see that overall the QSGW is in good agreement with ex- 
periment but the 'hybridized' QSGIT+LDA Hamiltonian 
is in even better agreement. For G and B we calculate 
sclf-consistently the self-energy and charge density under 
the corresponding deformation. We found that if instead 
we use the self-consistent self energy of the undistorted 
structure the value of G differs from that presented in Ta- 
ble I by «5%. In these calculations the atomic positions 
were allowed to relax within LDA in order to account 
for the displacement of the anion and cation sublattices 
relative to each other. For a pure shear deformation in 
the [111] direction this displacement can be viewed as 
a length change Al of the [111] bond described by the 
internal strain parameter £, Al = 3(1 — £)e/\/4- Our 
calculated value of £ is 0.53, in good agreement with pre- 
vious calculations 33 ' 34 

We apply two different deformations, the first of which 
is described by the following strain tensor: 

fe l e 0\ 
e = e ex (4) 
V e 2 J 

where e x =0.0025186, e 2 =-0.0049628 and e =0.0074814. 
This tensor conserves the volume and induces the B re- 
lated term in Vl l stress . To separate this term from the 
G related term we also performed a calculation with a 
deformation described by the following strain tensor: 

f° 6 °\ 

e = e (5) 

\ 0/ 

This strain tensor conserves the volume only to first or- 
der of deformation but only the first term (G related) is 
present in ni tress . 
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FIG. 2: (a) The shift of the conduction band minimum away 
from the T point in GaAs under deformation given by Eq. 
2] ('hybridized' method), (b) The magnitude of the conduc- 
tion band splitting along the [010] direction for the case of 
the 'hybridized' Hamiltonian, red line with deformation Q 
black line with deformation © (c) same as (b) but for QSG1T 
Hamiltonian. 

In Figs.^b) and (c) we show the magnitude of the split- 
ting along the [010] direction for the case of 'hybridized' 
and QSGW Hamiltonians, respectively. The black solid 
lines are for strain (0| and the red solid lines for ([5]) . The 
slope of the red line is equal to 

si = ^(Ce) 2 + B 2 (e 2 -ei) 2 (9) 
while the slope of the black line is equal to 



s 2 = |C| |e| (10) 

III. RESULTS 



A. Spin splittings 

In Fig. [2] we show the fc-dependence of the conduction 
band splitting along [010] for the case of 'hybridized' and 
QSGW calculations. Figure [DJa) shows the energy dis- 
persion around the T point for the case of the 'hybridized' 
Hamiltonian and the deformation given by Eq. 0J Along 
the [010] direction Sl^ vanishes and the dispersion is given 
by 

E±(k) = h 2 k 2 /2m eff ±1/2 \A\k (6) 

where 

\A\ = V(Ce) 2 +£? 2 (e 2 -ei) 2 (7) 

and the strain components were introduced in Sect. II. 
Correspondingly, the conduction band minimum shifts to 

fc± =±m eff /2h 2 \A\ (8) 



Hence 

\C\=s 2 /\e\ (11) 

and 

\B\ = y/al-4/\e 2 -e 1 \ (12) 

The magnitudes of C and B extracted with this proce- 
dure are given in Table I. As expected, in the case of 
the 'hybridized' method they are slightly larger than in 
the QSGW because the 'hybridized' band structure has 
a smaller band gap. However, the ratio of C / B remains 
nearly constant; it is equal to 3.197 in the 'hybridized' 
method and 3.171 in QSGIT. The magnitude of C is 
in good agreement with experiments in both the QSGW 
and the 'hybridized' methods but the latter should be 
trusted more due to better agreement of the other band 
parameters with experimental values. We also note that 
in the experimental determination of C, the B terms were 



4 



> 6 
D 

^4 
W 

^2 





^ 2 

> 
u 

1 



1 1 1 1 1 1 
- (a) e > 

\ \ — | — i — h- 


1 i > 
' 1 i 1 J," 


. (b) e < 




i.i.i 





0.02 0.04 0.06 0.08 
k[110] (A 1 ) 



0.1 



FIG. 3: (a) The conduction band splitting along the [110] di- 
rection for the case of 'hybridized' Hamiltonian, applied strain 
([5]) for e > (c) same as (a) but for e < 0. 



considered to be negligible. This corresponds to extract- 
ing the value of C directly from si; according to our cal- 
culations this would yield an inaccurate value of C = 7.14 
eV-A with an error of 5%, much less than the experimen- 
tal error in Ref^. The fl z D is highly anisotropic and 




for a general direction. It is therefore interesting to com- 
pare with Q l stress . We start by comparing C to 7. 
When we apply deformation (O, the dispersion along the 
[110] direction is 

E(k) = H 2 k 2 /2m eff ± 1/4 (7/fc 3 + Cek) (13) 
Thus the Dresselhaus and the stress terms can either add 
or subtract, depending on the relative sign of Ce and 
7. If they subtract the splitting will be zero at k rev = 
\J C/^y^/e i.e. the spin splitting reverses its sign at k rev . 
In Table I it is seen that y/Cpy = 0.895A -1 , which means 
that for the deformation ([5|) the spin splitting along [110] 
will change sign at k rev = 0.077 A~ x . As shown in Fig. 
1, e can be either negative or positive, therefore such 
cancellation will always occur depending on the sign of 
the C ■ e ■ k product. For example, if k > 0, whether 
it occurs for e > or e < depends on the sign of 
C. In Ref. [17J we determined the sign of 7 according to 
conventions in Ref. E Here we will determine the sign 
of C relative to the sign of 7 by simply plotting AE along 
[110] for positive and negative e. In Fig. [3] we show such 
plot for the deformation ([5]) with positive and negative 
e. The splitting is linear only in the vicinity of the V 
point, away from the T point the cubic term is clearly 
visible. For e > the two contributions add (Fig. [3ja)) 
but for e < they oppose each other (Fig. [3](b)); for 
k > 0.074A~ 1 the cubic term dominates and the splitting 
becomes positive. It is clear that C and 7 have the same 
sign (according to the convention used here they are both 
positive) . The sign of B is defined in a similar way. We 
apply the deformation: 



(14) 



So that £2 = — 2ei. Then the splitting along the [110] 
direction is AE oc [jk 3 + (ei — t2)Bk\ . if AE crosses 
zero for (ei — £2) > then B and 7 have opposite signs, 
otherwise B and 7 have the same sign. As can be seen 
in Fig. H] we find that B > 0. 
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B. Spin relaxation 



FIG. 4: (a) The conduction band splitting along the [110] 
direction for the case of 'hybridized' Hamiltonian, under de- 
formation (fl4)) with (£1 — 62) > (b) same as (a) but for 
(ei -e 2 ) < 0. 

completely vanishes in certain directions but is present 



After having reliably determined the values of material 
parameters that dictate the spin relaxation rate in GaAs 
it will be interesting to estimate the spin relaxation time 
for a deformation like the one given by Eq. [U In the 
Appendix we have derived the average spin relaxation 
time when both J~2d and f2 str ess are present 
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(Cef + (BAef 



J - 7l BCeAe-2-(v 2 T p (E)) (16) 



— ="fi 
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(17) 



Where || denotes axis parallel to the vector N = (0,0, 1) 
and _L perpendicular to it. Namely, (_L, 1) and (_L, 2) 
are the axes along the [110] and [110] crystal direc- 
tions, respectively. We see that there is an in-plane 
anisotropy induced by the simultaneous presence of e and 
Ae = ei — £2 strain components (a similar anisotropy 
was observed for the circular piezo-birefringence and 
confinement-induced circular birefringence in GaAs^). 
We can write the strain tensor (TH| in these axes: e = 
CijSij with ei,i = e[n ] =£i + £, e 2 ,2 = £[no] = £i - e and 
e 3,3 = e [ooi] = e 2- If we apply a uniaxial pressure (stress), 
p, along [110] in this system of coordinates, then using the 
compliance constants Sn, Si 2 and 644 we find £[110] — 
(S n + S12 + 5 44 /2) p/2, £[ii 0] = (S u + S 12 - S AA /2)p/2 
and £[001] = 2Si2p/2, hence, Ae = (Sn — S±2)p/2 and 
£ = (l/2)S 4 4P/2. According to Eqs. OS] and Q1 the dif- 
ference between the inverse of the spin relaxation time 
along [110] and [110] is equal to 

1 1 

T±,l T|_,2 

^H^(v 2 T p (E))Sn (Sn - S l2 )p 2 (18) 

Hence, an experimental setup similar to the one described 
above should be able to measure a linear increase of the 
difference (|18[) with the square of applied pressure. The 
rate of increase should be proportional to BC. In the 
experiment of Ref. 03 the authors measured the increase 
of the spin relaxation time with applied pressure along 
[100]. Such strain will only induce the B terms in Eqs. 
(fT5)) - (|17p . If we assume that the applied strain is large 
enough to ignore the Dresselhaus term, then the spin re- 
laxation time should be isotropic and should increase lin- 
early with the square of applied strain. However, unlike 
the experimental setup proposed here the rate of increase 
is proportional to B 2 . The experiment proposed here is 
independent of Dresselhaus terms no matter how small is 
the deformation, also the linear increase is proportional 
to CB instead of B 2 , hence it may be easier to detect. 
Provided that the orientation of the As-Ga bond has been 
previousely determined, this experiment can be used to 
find the sign of B relative to that of C from the sign of 
the difference (fI5|). 

IV. CONCLUSION 

We have presented first principles calculations of the 
magnitude and sign of bulk constants that govern the DP 



spin scattering in GaAs under strain. To our knowledge, 
this was the first estimation of magnitude and sign of B. 
We find that both C and B have the same sign as 7. Our 
value of C is in good agreement with experiments. We 
have derived an expression for the spin relaxation time 
of electrons under a strain given by Eq. U and showed 
that the in-plane spin relaxation is anisotropic in this 
case. We proposed an experiment that can exploit this 
anisotropy to deduce the magnitude and sign of B. 
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V. APPENDIX: SPIN SCATTERING RATE 

The momentum dependent spin relaxation time tensor 
is defined as: 

-L_ = llTp{E) (jW-^) (19) 

and 

-^^ llTp (E)(nf}-) (i^j) (20) 

Here i = x,y, z and the overbar denotes averaging over 
all directions of k. t p (E) is the momentum scattering 
time for an electron with energy E and 

f^a(cosd)(l - Pi (cos 9) )d cos 9 

7i = +1 (21) 

J , a (cos 9) (1 ~ cos 9) d cos 9 

where a(cos9) is the electron scattering cross section and 
Pi the Legendre polynomials. Here it is assumed that the 
electron scattering is elastic, the electron energy spec- 
trum is isotropic and the scattering cross section er(k, k') 
depends only on the scattering angle 9. ft is the total 
effective field and in our case we can write: 

ft = flrj + J~2 stress (22) 

with the components as given in Eqs. fl} and @. We 
apply a strain like that of tensor Eq. (H} with the con- 
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straint e 2 = — 2ei, so as to conserve the volume. Then 
^stress = Ce^y + BAek x 

[Cek x + BAek y ] (23) 



stress 



^stress ^ 



where Ae = (ei — £2)- To facilitate the discussion let's 
write also explicitly the components of the Dresselhaus 
field: 



QX 


= 2-/k x (fc 2 


-kl) 




li D 


= Zlky {kl 


-kl) 


(24) 


O 2 


= 2 7 fc z (k 2 x 


~kl) 





Then we can write 



^x — (^stress) 2 + (^L>) 2 + ^stress^D 

The first integral on the RHS is 



(25) 



(ft* tress ) 2 = {Ce ) 2 kl + (BAe) kl + (2CBeAe) k x k v 



(Ce) + (BAe) 



The second integral on the RHS is 



16 



(^)2= 4 7 2 fc2(fc2_fc2)^ =7 2fc6 



105 



(26) 



(27) 



The third integral on the RHS is 



^stress^D = (^l) kxky (k 2 y - fc 2 ) + 

+ (BAe 1 )k 2 x (k*-k*)=0 



So we get 



3 



(Ce) 2 + (BAe) 2 



fc 2 + 7 2 fc 6 i| (28) 



In a similar way we obtain £l y = Q x and 



105 



For the off-diagonal components we get 



(29) 



OQ — OO — Q x Q y 4- Q x O y 4- 
I 755 7yy~ , 757^ Tyy 

^ iL stress iL D ~ J ' stress* ^stress 



(30) 



The first three integrals on the RHS are equal to zero. 
With the form of strain field given by Eq. (23) the last 
term can be written as 



QX QV 

iL stress iL stress 



-(Ce) 2 k y k x + (BAe) k x k y - CBeAe(k 2 + k 2 y ) = 



-CBeAek 2 



(31) 



All other f^fl,- are equal to zero. Therefore according to 
equations (1) and (2) for the spin relaxation time of an 
electron with energy E we find 
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- lir P (E) 
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(32) 
(33) 
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Then the average spin relaxation time is 
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32 2 TO 6 
105 7 F 
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(Cef + (BAef 



(v\(E)) 



32 9 m 6 
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(36) 



1 1 2 m 2 

= = -~1iBCeAe-^{v\ p (E)) 



(37) 



where v = Tik/m and the brackets ( ) denote averaging By transforming the above tensor to the principal sys- 
over energies. For example for the Maxwell distribution tern of coordinates we obtain: 



— = ^111 (Cef + (BAef 



32 o to 



' <Ap(£)) + i05 7 V ( 



rr> 



+ -7 l BCeAe- s -{v*T J> (E)) (38) 



Tl 2 



= H 



[Cl r + iBA, r ^(v\(E)) + ^-^( V \(E)) 



- -jtBCeAe^-iv^E)) (39) 
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r ll 



2 r 

3 
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32 9 to 



{Cef + (BAeY -^(v\(E)) + _y _(,,%(£)} 



(40) 



Where || denotes axis parallel to the vector N = (0,0, 1) 
and _L perpendicular to it. Namely, (_L, 1) and (_L, 2) 
are the axes along the [110] and [110] crystal directions, 



respectively. Equations (f3"8"]l and (|3"9"|) signal an in-plane 
anisotropy induced by the simultaneous presence of e and 
Ae. 
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